DESY 95-210 
KA-TP-8-95 
hep-ph/9511344 
November 1995 



QCD Corrections to Hadronic Higgs Decays 

A. Djouadi 1,2 Q, M. Spira 3 and P.M. Zerwas 2 

1 Institut fur Theoretische Physik, Universitat Karlsruhe, 
D-76128 Karlsruhe, FRG. 

2 Deutsches Elektronen-Synchrotron DESY, D-22603 Hamburg, FRG. 

3 II. Institut fur Theoretische PhysikQ D-22761 Hamburg, FRG. 



Abstract 

We present an update of the branching ratios for Higgs decays in the Standard 
Model and the Minimal Supersymmetric extension of the Standard Model. 
In particular, the decays of the Higgs particles to quark and gluon jets are 
analyzed and the spread in the theoretical predictions due to uncertainties of 
the quark masses and the QCD coupling is discussed. 
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1. Introduction 



The coupling of the Higgs bosons to other particles grows with the mass of the particles. 
This characteristic property is a direct consequence of mass generation through the Higgs 
mechanism. To establish the Higgs mechanism experimentally, it is therefore mandatory 
to measure the couplings very accurately Q once scalar particles have been found. The 
main test grounds for the Higgs couplings to gauge bosons are the production cross sections 
for Higgs-strahlung off gauge bosons and WW/ZZ fusion, and the widths/branching ratios 
for Higgs decays to gauge bosons. The Higgs couplings to heavy quarks determine the 
cross sections for the production of Higgs particles in gg fusion at hadron colliders [§, [| , 
as well as the rate of Higgs bremsstrahlung off heavy quarks at e + e _ |3J] and hadron 
colliders || . The measurement of Higgs decay branching ratios, including b, c quarks and 
t leptons 0, provides a complementary method to determine the Higgs couplings. 

In this note we will reanalyze the branching ratios for Higgs decays to b, c quark 
jets and to light hadron jets evolving out of gluon decays, 

H -> bb/cc + ... (1) 
H - gg+ ... (2) 

The ellipses indicate additional gluon and quark partons in the final state due to QCD 
radiative corrections. Special attention will be paid to uncertainties related to the b, c 
quark masses and the QCD coupling a s . It turns out that the evolution of the charm 
quark mass from low energy scales, where it can be determined by QCD sum rules, to 
high energy scales defined by the Higgs mass, introduces very large uncertainties in the 
cc branching ratio. The partial width of the second decay mode (0) will be derived for 
gluon and light quark final states since heavy quarks add to the partial width of the 
first decay process ([!]). The b, c and gluon decay modes are experimentally important in 
the Standard Model (SAi) for Higgs masses less than about 150 GeV. In the minimal 
supersymmetric extension (M.SSM) b quark decays may be dominant for a much wider 
range in the parameter space. 



2. Standard Model 

2.1 b, c quark decays of the SM. Higgs particle 

The particle width for decays to (massless) b, c quarks directly coupled to the SM. Higgs 
particle is given, up to 0(a 2 s ) QCD radiative corrections^], (Fig. la) by the well-known 
expression M, 



10 



Of JUT 

T[H - QQ] = -^=± m 2 Q (M H ) [A QCD + A t ] (3) 



1 Thc effect of the efectroweak radiative corrections in the branching ratios is negligible [fTD. 
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in the MS renormalization scheme; the running quark mass and the QCD coupling are 
defined at the scale of the Higgs mass, absorbing this way any large logarithms. The quark 
masses can be neglected in general except for top quark decays where this approximation 
holds only sufficiently far above threshold; the QCD correction in this case are given in 
the Appendix. 

Since the relation between the pole mass M c of the charm quark and the MS mass 
evaluated at the pole mass m c (M c ) is badly convergent ]T2|, we will adopt the running 
quark masses WIq(Mq) as starting points. They have been extracted directly from QCD 
sum rules evaluated in a consistent 0(a s ) expansion [|H|]. The evolution from Mq upwards 
to a renormalization scale /i is given by 



m Q 0) = m Q (M Q ) 
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for M c < [i < M b 
for M b < fi 



For the charm quark mass the evolution is determined by eq.(Q) up to the scale \x = Mb, 
while for scales above the bottom mass the evolution must be restarted at Mq = Mb. 
Typical values of the running b, c masses at the scale /i = 100 GeV, characteristic for the 
Higgs mass, are displayed in Table 1. The evolution has been calculated for the QCD 
coupling 



a s {M z 



0.118 ± 0.006 
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defined at the Z mass ||14|| . The large uncertainty in the running charm mass is a conse- 
quence of the small scale at which the evolution starts and where the errors of the QCD 
coupling are very large. 



2.2 Higgs decay to light hadron jets 

The decay of the Higgs boson to gluons is mediated by heavy quark loops in the Standard 
Model (Fig. lb); the partial decay width JOJ is given by 



T LO [H -> gg] 
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a s {M z ) 


m Q (M Q ) M Q = M% 2 


m Q (fj, = 100 GeV) 


b 


0.112 
0.118 
0.124 


(4.26 ± 0.02) GeV (4.62 ± 0.02) GeV 
(4.23 ± 0.02) GeV (4.62 ± 0.02) GeV 
(4.19 ± 0.02) GeV (4.62 ± 0.02) GeV 


(3.04 ±0.02) GeV 
(2.92 ±0.02) GeV 
(2.80 ± 0.02) GeV 


c 


0.112 
0.118 
0.124 


(1.25 ± 0.03) GeV (1.42 ± 0.03) GeV 
(1.23 ± 0.03) GeV (1.42 ± 0.03) GeV 
(1.19 ± 0.03) GeV (1.42 ± 0.03) GeV 


(0.69 ± 0.02) GeV 
(0.62 ±0.02) GeV 
(0.53 ± 0.02) GeV 



Table 1: The running b, c quark masses in the MS renormalization scheme at the scale 
\x = 100 GeV. The starting points tTIqIMq) of the evolution are extracted from QCD 
sum rules J73|// the pole masses Mq 2 are defined by the 0(a s ) relation mQ(MQ 2 ) = 
Mq 2 /[I ± Aa s /37f] with the running masses. 
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The parameter tq = AMq/Mfj is defined by the pole mass of the heavy loop quark Q. 
For large quark masses the form factor approaches unity. QCD radiative corrections are 
built up by the exchange of virtual gluons and the splitting of a gluon into two gluons 
or a quark-antiquark pair, Fig. lb. If all quarks u,---,b are treated massless at the 
renormalization scale fi ~ Mu ~ 100 GeV, the radiative corrections can be approximated 
very well by (16 



r"' [H - 99(g), qqg] = r LO [a<"'>00] ii + E"' 
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with Np = 5 light quark flavors. The radiative corrections are very large, nearly doubling 
the partial width. 
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The final states H — ► bbg and ccg are also generated through processes in which the 
b, c quarks are coupled to the Higgs boson directly (Fig. la). Gluon splitting g — » bb in 
H — > gg (Fig. lb) increases the inclusive decay probabilities0r(/7 — > 66+. . .) etc. Since b 
quarks, and eventually c quarks, can in principle be tagged experimentally, it is physically 
meaningful to consider the particle width of Higgs decays to gluon and light u, d, s quark 
final jets separately. The contribution of b, c quark final states to the coefficient E F in 
eq.(0) is given by 
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Instead of naively subtracting this contribution, it may be 



in the limit M\ > 

noticed that the mass logarithms can be absorbed by changing the number of active 
flavors from N F = 5 to N L = 3 in the QCD coupling, 
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This way we arrive again at an equally simple expression 



r[/f-> gg(g),gqg] = T LO 
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for Nl = 3 light q = u,d,s quark flavors in the final state. 

The subtracted parts may be added to the partial decay widths into c and b quarks. 
Up to 0(a 3 B ), they are given by the difference of the gluonic widths [eq.(|7])] for the corre- 
sponding number of flavors iV^, 
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The MS A parameters for three, four and five flavors in the QCD coupling 
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(10) 



are given in Table 2 together with the effective couplings a^ Np '(Mz) for three, four and 
five flavors corresponding to the QCD scales A]y F . In a^(Mz) the contribution of the 6 
quark is subtracted and in af'(M^) the contributions of both the b and c quarks. 



2 The two contributions add up incoherently in the limit where the final state quark masses are ne- 
glected [apart from the Yukawa Higgs couplings] . The topology of the final states is in general different 
for decays to b, c quarks through the direct couplings or gluon splitting: in the former decay mode, quark 
and anti-quark jets are emitted primarily back-to-back, while in the latter decay mode they tend to be 
collinear. 



5 



af\M z ) 


A 5 [MeV] 


«i 4) (M z ) 


A 4 [MeV] 


af\M z ) 


A 3 [MeV] 


0.112 


159 


0.107 


238 


0.101 


286 


0.118 


226 


0.113 


327 


0.105 


378 


0.124 


312 


0.118 


434 


0.110 


483 



Table 2: QCD scales An f for the range of uncertainty in the strong coupling constant 
otf\Mz) = 0.118 ± 0.006. The effective couplings a^'(Mz), with the b quark decoupled, 
and af\Mz), with both the b and c quarks decoupled, are given for comparison, too. 



With E 3 = 20.25, the QCD radiative corrections still amountQto ~ 70%. However, a 
large fraction of the corrections can be absorbed by choosing, for the proper renormaliza- 
tion scale, the BLM scale |T7| which maps contributions associated with gluon self-energies 
into the effective QCD coupling; this is technically implemented by choosing /x such that 
the coefficient of the Np or jVjr term vanishes: 
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The QCD corrections to the partial width 

r[# -> gg(g),qqg] = T LO 
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are reduced in this approach to a comfortable level of 15 to 25%. 



2.3 Numerical evaluation 

The numerical analysis of the branching ratios for the Higgs decays in the Standard Model 
has been performed for the set of parameters given in the tables and the top quark massQ 

M t = (176 ± 11) GeV 

To estimate systematic uncertainties, the variation of the c mass has been stretched over 
2<j and the uncertainty of the b mass to 0.05 GeV. However, the dominant error in the 
predictions is due to the uncertainty in a s which migrates to the running quark masses 
at the high energy scales. 

The results for the branching ratios are displayed in Fig. 2. Separately shown are 
the branching ratios for r's, c, b quarks, gluons plus light quarks and electroweak gauge 

If the b, c quarks are included in the final state, the partial width increases by an additional 20%. 
4 This value of the top quark mass M t corresponds to the average of all measurements at the Tevatron 
presented in Ref . p8[ . 
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bosons. The uncertainties in the prediction for the charm and gluon branching ratios are 
very large. Increasing a s reduces the value of the running c mass quite dramatically^. 

Nevertheless, the expected hierarchy of the Higgs decay modes is clearly visible in 
Fig. 2 despite these uncertainties. BR T is more than an order of magnitude smaller than 
BRb, a straight consequence of the ratio between the two masses squared and the color 
factor. As a result of the small charm quark mass at the scale of the Higgs mass, the ratio 
of BR C to BRb is reduced by much more than an order of magnitude, which would have 
been naively expected. Taking these subtle QCD effects into account, the measurement 
of the decay branching ratios provides an excellent method to explore the physical nature 
of the Higgs particle. 



3. The Minimal Supersymmetric Standard Model 



We have performed a similar analysis for the hadronic decay modes of the Higgs bosons 
h, H, A, H ± in the Minimal Supersymmetric extension of the Standard Model (AiSSAi). 
Apart from the usual modifications Qq of the couplings, the analytic expressions for the 
partial widths of the scalar neutral Higgs bosons h, H are the same as in the Standard 
Model, eqs.(^J) and (0). In the massless quark limit [the general massive case is treated 
in the Appendix], the QCD radiatively corrected decay widths into quarks are given by 
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with $ = h,H, A. [Eq. (p!3|) is valid if either the first or the second term is dominant.] The 
relative couplings gq have recently been collected in Ref.|20]; the masses in the Yukawa 
couplings are to be evaluated at the scales M$ and Mjj±- 

Since the b quark couplings to the Higgs bosons may be strongly enhanced and the t 
quark couplings suppressed in the M. SSM., b loops can contribute significantly to the gg 

5 The value of BR C is significantly smaller in Ref.0 for two reasons, (i) The ratio between the MS mass 
m c (M c ) and the pole mass M^' 3 to 0{ct 2 s ) is smaller than the corresponding ratio for M^* 2 to 0(ce s ). 
However, since m c (M c ) has been determined through QCD sum rules only to 0(a s ), the parameter 
M c = MP t2 is the proper pole mass to be used in a consistent analysis up to 0(a s ). We have performed 
the evolution of the running MS mass with and without the 0(a 2 s ) contribution; the difference between 
the two results at the scale of the Higgs mass turned out to be negligible. The present analysis is therefore 
theoretically consistent, (ii) Moreover, in Ref. the average LEP a s value has been adopted which is 
larger than the world average value including deep-inelastic scattering data. This gives rise to a faster 
fall-off of the running charm mass ~ [a s (^)] 12 / 23 at large scales. 
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coupling so that the approximation Mq 3> Mfj cannot be applied any more in general. 
Nevertheless, it turns out a posteriori that this is an excellent approximation for the 
QCD corrections. The LO width for h, H — > gg is obtained from eq.(|6|) by substituting 
A H — > gg 11 A h,H ; for the pseudoscalar Higgs decays, we find H 



T LO [A -> gg] 
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T[A ^ gg(g),qqg] 
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To illustrate the size of the uncertainties introduced into the predictions by the QCD 
parameters, the branching ratios have been calculated for a specific set of parameters. 
The top mass is varied within M t = (176 ± 11) GeV. In addition to the other parameters 
defined in the previous section, the running mass of the s quark at the scale 1 GeV and 
the CKM type mixing parameter V c b are chosen as 

m s (l GeV) = (0.190 ± 0.040)GeV 
\V cb \ = 0.040 ±0.008 

while the SUSy parameters are set to 

tan/3 = 1.6 

SUSy scale parameter: M s = 1 TeV 
"typical" mixing: A t = —fx = 1 TeV 

SUSy masses and couplings have been calculated according to the RG program described 
Varying the SUSy parameters does not change the picture of the QCD 



in Ref. 19 



corrections and the uncertainties exemplified for the set of parameters chosen above. The 
final results^ are displayed in Figs.3a-d. The branching ratios are separated again for 
final states including b, c quarks [labeled bb and cc] and gluons plus light quarks [labeled 
gg]. The main sources of uncertainties for the branching ratios are the charm and gluon 
decays. The branching ratios for b and r decays are less affected by a s . The uncertainty 
in the top quark mass affects primarily the upper limit of the light CP-even scalar mass 
Mh and the couplings at the electroweak level, as shown in Fig.3a. The uncertainty in 
M h migrates to the partial widths of the heavy Higgs bosons through cascade decays; this 



6 The Fortran code for the partial decay widths in the SM. and MSSM may be obtained from 
djouadi@desy.de or spira@desy.de. 
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second step is indicated by the error bars attached to the curves in Figs. 3b-d. 
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APPENDIX 



For completeness, we present in this Appendix the expressions of the leading order QCD 
corrections to SAi and A4SSA4 Higgs boson decays involving non-zero mass effects for 
heavy quarks. As a general example we will consider top quarks. 

The partial decay widths of the CP-even Higgs bosons $ = Hsm , h and H into top 
quark pairs, in terms of the top quark pole mass, is given by 



r[$ -> tt] 
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where j3 — (1 — 4m^/M|) 1/ ' 2 is the velocity of the top quarks. To leading order, the QCD 
correction factor is given by H, 121 
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[Li 2 is the Spence function defined by Li 2 (x) = — Jjf rfyy 1 log(l — y).] 

The partial decay width of the CP-odd Higgs boson A into top quark pairs reads 
correspondingly 



T\A^tt] = Z -^ {!l tfmlf> 
where the QCD correction factor is given by [JO, E!3 
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The partial decay width of the charged Higgs particles decay into top and bottom 
quarks may be written 
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where //£ = mj/M^, fib = ml/M^j± and A the usual two-body phase space function A = 
(1 — fit — fib) 2 — 4fi t fib] again the quark masses are the pole masses. In the approximation 
where the b quark mass is neglected the QCD factors A^ are given by |22|] 
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Well above the tt threshold, the QCD factors Aj reduce to 



A 



tb 



A 



bt 



9 3, 

- + - loe 

4 2 £ 

9 3, 

- + - loe 

4 2 ^ 



m t 



(A.7) 



The large logarithms can be mapped into the running quark masses in the usual way. 
Adopting the MS mass at the scale of the Higgs mass the bulk of the next-to-leading 
order correction is automatically included in this limit and the QCD corrections approach 
the common chirally invariant factor Aqcd, ec i-(@)- 

In the small mass limit, the non-logarithmic term of the factor A^ 
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differs from the corresponding term in A = t,b), eq.( |A.7[) . However, this is still in 

accordance with chiral symmetry since the correction eq.( [A.8|) is of subleading order in 
the small quark mass expansion. 
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Figure Captions 

Fig. 1: (a) Generic Feynman diagrams for the decay processes H — > bb,cc and H — > 
bbg, ccg. (b) Feynman diagrams for H —> gg decays and final-state gluon splitting into 
quarks. 

Fig. 2: Branching ratios of the SAi Higgs boson including the uncertainties from the 
quark masses and the QCD coupling a s . The resulting errors of the branching ratios are 
presented as shaded bands. The curves labeled bb and cc are the inclusive decay widths; 
they account for all final states including c and b quarks. The curve labeled gg corresponds 
accordingly to gluon and light-quark final states only. 

Fig. 3: Branching ratios of the M.SSM. Higgs bosons h, H, A, H ± [(a). .(d)], including 
uncertainties from the quark masses mj, m c , m s and the strong coupling a s . In the widths, 
the three-body channels [ED| have been included. The top quark mass is fixed at M t = 



176 GeV for the shaded error bands and the error bars shown below the 200 GeV mark in 
figure (b). The additional uncertainty due to the top mass is marked by error bars in the 
figures (b) [at M H > 200 GeV], (c) and (d). In figure (a) the curves for the upper and 
lower limit of the top mass band are presented separately, using the average values of the 
other quark masses and of the strong coupling a s for the sake of clarity. The labels follow 
the definitions in Fig.2; i.e. the branching ratios are classified according to the inclusive 
hadronic final states with [labels bb, cc] and without heavy quarks [label gg] . 
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